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Abstract 

We show that many functions containing W are Stieltjes functions. Explicit Stieltjes integrals are 
given for functions 1/W{z), W{z)/z, and others. We also prove a generalization of a conjecture of 
Jackson, Procacci & Sokal. Integral representations of W and related functions are also given which 
are associated with the properties of their being Pick or Bernstein functions. Representations based 
on Poisson and Burniston-Siewert integrals are given as well. 

1 Introduction 

The Lambert W function is the multivalued inverse of the mapping W t-^ We^ . The branches, denoted 
by Wfc (fc G Z), are defined through the equations [I] 

VzeC, VFfe(z)exp(l¥fe(z)) = 2, (1) 
VFfe(z) ^> lufc z as — > oo , (2) 

where In/j z = \nz + 27rifc, and Inz is the principal branch of natural logarithm [2 . This paper considers 
only the principal branch k — 0, which is the branch that maps the positive real axis onto itself, and 
therefore we shall usually abbreviate Wq as W herein. 

The W function is rich in integral representations. As a consequence, is a member of a number 
of function classes, specifically, the classes of Stieltjes functions. Pick functions and Bernstein functions, 
including subclasses Thorin-Bernstein functions and complete Bernstein functions. A description of the 
classes can be found in a review paper [3] and a recently published book [1] . The integral representations 
are related to the fact that is a real symmetric function, in the terminology of [3 p. 160] (see also [SI 
p. 155]), with positive values on the positive real line. In addition to integral representations connected 
to the above function classes, integrals following Poisson [33] and Burniston-Siewert [5] are given. 

The classes of Stieltjes functions and Bernstein functions are tightly connected with the class of 
completely monotonic functions that have many applications in different fields of science; the list of 
appropriate references is given in [3]. Therefore we shall also study the complete monotonicity of some 
functions containing W. 

The properties and integral representations mentioned above have interesting computational implica- 
tions. For example, that W{z)/z is a Stieltjes function means that the poles of successive Pade approx- 
imants interlace and all lie on the negative real axis [3 p. 186] (here in the interval — oo < z < — 1/e). 
In addition, some of the integral representations permit spectrally convergent quadratures for numerical 
evaluation. 



1.1 Properties of W 

For convenience, we recall from W some properties of W that are used below. The function is discontin- 
uous on the branch cut B C K, defined to be the interval B = (— oo, —1/e]. On the cut plane C\B, the 
function is holomorphic. Its real values obey —1 < W{x) < for a; G [—1/e, 0), W{0) — and W{x) > 
for a; > 0. The imaginary part of W{t) has the following range of values for real t. 

^W{t) e (0, tt) for t e (-00, -1/e) and ^W{t) = otherwise. (3) 
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QW{z) and have the same sign in the cut plane C\M, or equivalently 

3H/(z)3z > . (4) 



W has near conjugate symmetry, meaning W{z) — W{z) except on the branch cut B. We also note 

Wiz) =\nz-\nW{z) (5) 
in the cut plane C\(— oo, 0]. The Taylor series near z = is 

Wiz) = J2{~nr~''- (6) 
with radius of convergence 1/e, while the asymptotic behaviour of W{z) near its branch point is given by 



W{z) - -1 + v/2(ez+ 1) as z ^ -1/e . (7) 
It follows from ^ and © that 

W{z)/z^l as z^O. (8) 
W{z)/z ^ as z oo . (9) 

li z = t + is and W{z) = u + iv, then 

e"(ucosw — usinw) — t, e^{usmv + vcosv) — s . 

For the case of real z, i.e. s = 0, the functions u — u{t) and v — v{t) are defined by 

u — —V cot V, (10) 
t = t{v) = -i;csc(i;)e-"'=°*'' . (11) 

For the case of purely imaginary z, i.e. t = 0, the functions u = u{s) and v = v{s) obey 

u vtaiiv, (12) 
s = s(i;) = t-sec(t;)e^'^"V (13) 



z(l + W{z)) 



W'{z) ^1 as z ^ , (15) 
VF'(z) ^0 as z oo . (16) 



dW(z)/dz^dW{z)/dz (17) 



The derivative of W{z) is given by 
It follows from © and 

Near conjugate symmetry implies 
for z ^ (-0O, 0]. 

1.2 Stieltjes functions 

We now review the properties of Stieltjes functions, again concentrating on results that will be used in 
this paper. We must note at once that there exist several different definitions of Stieltjes functions in the 
literature, and here we follow the definition of Berg [3]. 

Definition 1.1. A function f : (0, oo) — > M is called a Stieltjes function if it admits a representation 

fix)^a+ for x>0, (18) 

Jo x + t 

where a is a non-negative constant and fi is a positive measure on [0, oo) such that (l + t)^^ dij.{t) < oo. 
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A Stieltjes function is also called a Stieltjes transform [9, p. 127], and the integral in the right-hand 
side of is called a Stieltjes integral of f{x). Except in H2.'3\ below, the term Stieltjes function will 
here always refer to Definition [TBI 

Theorem 1.2. The set S of all Stieltjes functions forms a convex cone fW, p. 127] and possesses the 
following properties. 

(i) / e5\{0}^^e5 

(u) /e5\{0}^^e5 

(tit) f^S^^^S (c > 0) 

(iv) /,5e5\{0}^/oi e5 

(v) /,5e5\{0}^^e5 

(vi) f,geS ^ /".g^"" (0 < a < 1) 
(vii) / e 5 ^ /" e 5 (0 < a < 1) 



(vm) /e5\{0}=^i(S 



1] eS 



(ix) / e 5 \ {0} , linij;^o xf{x) = c > f{x) - c/x £ S 
(x) f eS^ f^{0)~ f"{l/x) eS (0<a<l) 

(xi) feS\{0}^^[l-j^)eS 

(xii) feS, lim.^oo fix) =c>0^{cP - f^)eS (-1 < /? < 0) 

Proof. Properties H]) are listed in [3]; property ([6]) is due to the fact that the Stieltjes cone is 

logarithmically convex [IV and property ([7]) is its immediate consequence. Property ([5]) is taken from 
[TUl p. 406]. Property (O follows from properties © and ([5]) in the following way: / G iS \ {0} ^ 
5(a;) = l/(a:/(a;)) € S ^ {giO)/g{x) - l)/x = ixf{x)/c ~l)/x e S ^ f{x) - c/x £ S. The last three 
properties PTI|) -(|12p will be proved in § [S] 



A Stieltjes function / has a holoniorphic extension to the cut plane C\(— oo, 0] satisfying /(z) = f{z) 
(see [n], and [1 p. 11-12]) 

f{z)^a+ (|argz|<7r) . (19) 

Jo z + t 

In addition, a Stieltjes function /(z) in the cut plane C\(— cxd, 0] can be represented in the integral form 

/W- / T-^ (|argz|<7r), (20) 
Jo i- + uz 

with a bounded and non-decreasing function $(«). The integral ([20| is useful for the analysis of Pade 
approximations of Stieltjes functions ^7, Ch. 5] and equivalent to the representation (I19|) . which follows 
from the following observation. 

According to properties Q and ([2]) , if a function / G 5 then f{l/x)/x G S as well and hence the 
latter admits representation ([TBI) 



which after replacmg x with 1/x gives 

dp{t) 



1 + a;i ' 



where the constant a can be included in the integral. Finally, one considers the holomorphic extension 
of the last integral to the cut plane C\(— oo, 0] similar to obtaining p^ . 

There are various kinds of necessary and sufficient conditions implying that a function / is a Stieltjes 
function. Some of them are based on the classical results established by R. Nevanlinna, F. Riesz, and 
Herglotz. Here we quote two such theorems, taken from [13l p. 93] and [3] Theorem 3.2]. 



Theorem 1.3. A function g{z) admits an integral representation in the upper half-plane in the form 

g{z) = I ^ (3- > 0) , (21) 

with a non- decreasing function $(w) of hounded variation on R (i.e. d'^{u) < oo for smooth ), if 
and only if g{z) is holomorphic in the upper half-plane and 

^giz) > and sup \yg{iy)\ < oo . (22) 

l<y<oo 

To apply Theorem [Ql to the integral ^ one should set g{z) = -/(-l/z)/z (cf. fT", (6.12) on p. 215]), 
then conditions ([22|) read as 

3/(-l/z)/z<0 and sup \f{i/y)\ <oo. (23) 

l<y<oo 

Theorem 1.4. A function f : (0, oo) — > R is a Stieltjes function if and only if f{x) > for a; > and 
there is a holomorphic extension f{z), z — x -\-iy, to the cut plane C\(— cxd,0] satisfying 

3/(^) < for 3z > 0. (24) 

Remark 1.5. The inequalities (|24p alone express a necessary condition for f to be a Stieltjes function. 
In the terminology of 110, p. 358], a holomorphic function f{z) is called a Herglotz function if 3/ > 
when 3z > 0, 3/ — when 3z = and SJ/ < when 3z < 0. Thus, for f to he a Stieltjes function it is 
necessary that f he an anti-Herglotz function (cf. ]1(A p. 406]). 



2 Stieltjes functions containing W{z) 

In this section we consider a number of functions containing W{z) and prove that they are Stieltjes 
functions. We begin with the function W{z)l z. 

2.1 The function W{z)/z 

Although we could establish the result conveniently by applying one of the criteria stated in > jl.21 we 
nonetheless first present a direct proof that is of great importance for further investigations. Moreover, 
compared with using the criteria above, the present way allows us to make useful observations which are 
given in the remarks following the proof and used in further evidence. 

Theorem 2.1. The function W{z)/z can be represented in the form of a Stieltjes integral (|20p . 

Proof. From ([8|), the function W{z)/z is single- valued and holomorphic in the same domain as W{z), 
namely D = {z £ C | z ^ B}, and can be represented by the Cauchy's integral formula 

- ' f -dt , (25) 



z 27rj Jq t{t — z) 

where C is the standard 'keyhole' contour which consists of a small circle around the branch point i = — 1/e 
of radius, say r, and a large circle around the origin of radius, say i?; the circles being connected through 
the upper and lower edges of the cut along the negative real axis. Then for sufficiently small r and large 
R the interior of the contour C encloses any point in D. 

Let us consider the integral ([25]) in the limit in which r -> and oo. Using asymptotic estimations 
([7|) and ([2]), it is easily seen that the contributions of each circle to the integral ([25]) go to zero. As a 
result, in accordance with the assignment of values of W function on the branch cut, the integral becomes 

^(z)_ 1 Wit) 



z 2'Ki J_rx, t{t — z) 27ri J^i/e t{t — z) 

which reduces to 



, t{t~z) 
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The change of variable t = —1/u transforms the integral to the form (PU)) 

w{z) _ r d<^iu) 

Z Jo 1 + uz ' 



(27) 



where 

= - / '^W{-l/t) dt, u e [0,e] . (28) 
7^ Jo 

According to (jS)) 

< 3W^(-l/t) < TT for i G [0,e] , (29) 
hence $'(m) > 0, and thus $(u) is a bounded and non-decreasing function. 

Remark 2.2. The Junction W{z)/z is a real symmetric function as is any Stieltjes function (this im- 
mediately follows from Definition 2.2), which just corresponds to the near conjugate symmetry property. 

Remark 2.3. It follows from the proof of the theorem that not only is W{z)/z a Stieltjes function, it 
has a smooth measure ^(u) defined by (|28p. Moreover, all the moment integrals 

f t"d'S>it) (n = 0,1,2,...) 
Jo 

exist because so do the integrals 

f"iW{~l/t)dt 







owing to (|29l) . This remark is useful for justifying the use of Fade approximants for the evaluation of 
W{z) mentioned in ^Tj. 

Now we apply Theorem ll.3l to the function f{z) = W{z)/z for which conditions (|23p read as 

3W^(-l/z) > and sup \yW{i/y)\ < oo . 

l<y<C30 

The first condition is satisfied by (|H) because 3(— 1/z) and 3z are of the same sign. To verify the second 
condition we set W{i/y) — u + iv and put s = 1/y in (|T2|) and ([T3)) . As a result, since < < 7r/2 for 
2/ > 0, we obtain 

\yW{i/y)f ^ y^u^ + v^) ^ y^v^il + Un' w) = y2„2/ ^032 ^ ^ g-2.tan„ < ^ 

To extend the result to the lower half-plane "^z < it is enough to take the complex conjugate of both 
sides of the representation ([20)1 and use the near conjugate symmetry of W. Thus Theorem [T3] gives us 
one more way to prove that W{z)/z is a Stieltjes function. 

2.2 Other functions 

By Theorem 12. 1[ W{x)/x € S. Using this result and the properties of the set S listed in section [L2l we 
now give some classes of functions that are members of S. 

Theorem 2.4. The following functions belong to the set S, for x > 0. 

(a) l/{c + Wix)), c>0 

(b) W°'{l/x),0<a<l 

(c) xl^W/^il/x),-! < 13 <0 

(d) W{x)/[xic + W{x))], c>0 

(e) l/Wix)-l/x 

(f) c+W{x^), c> 0, -1 < /3 < 

(g) l/{c+W{x°')),c>0,0<a<l 



(h) x°'^^W-"^{x^)[l + Wix^^)]^-^, 0<a<l, -l</3<0, 0<7<1 

(i) l/(c + a;"), < a < 1, c > 
(j) 1 - x°'W'^{l/x), < a < 1 

(k) 1 - x-°'^W°'{x'^)[l + Wix^)]-", 0<a<l, -1<^<0 
Proof. We use the properties listed in Theorem 11.21 

(a) We apply property ([2]) to H/"(a;)/x to find that l/VK(a;) e 5 and then apply dS]) to 1/W^(x). 

(b) We first apply ([U to f{x) — 1/W{x) that is in S by statement (a) and find W{l/x) £ S. Then we 
apply © to W{l/x). 

(c) Apply ([1]) to W{x)/x and apply then ([7]) to the result. 

(d) Apply (fTTj) to the function in the statement (a) using W{0) = 0. 

(e) Apply ^ to W{x)/x using ^ or apply (|9]) to the function in the statement (a) with c — 0. 

(f) Apply ^ to the function in the statement (a) and g{x) — x^ (—1 < /3 < 0) that is in iS [9l[3]. 

(g) Apply (|4]) to the function in the statement (a) and g{x) — x^" e for < a < 1. 

(h) Apply §1 to functions f{x) = W{x)/x and g{x) = a;^ (-1 < /? < 0) and find xl^W'^x^^) G S. 
Hence by © a{x) = x°''^W~°'{x'^) G 5 for < a < 1. Then apply © to the function in the 
statement (a) with c = 1 and g{x) ~ x^ to get b{x) = 1 + W{x^) G S. Finally apply ^ to a{x) 
and b{x). 

(i) It suffices to take © with c = 1 and f{x) = x^"' (0 < a < 1). 

(j) Apply to the function in the statement (c) with f3 — —1 using ([5]) (or apply (x) to W{x)/x). 

(k) Apply ([TU|) (or (IT^ ) to the result of application of (|1]) (respectively (0 ) to the function in the 
statement (d) with c = 1 and g{x) — x^ (—1 < /3 < 0). 

The next theorem proves and generalizes a conjecture in |17) . 
Theorem 2.5. T/ie following functions are Stieltjes functions for fixed real a 6 (0, e] 

Fo{z) = -^W{a{l + z))/[Wiail + z))-W{a)f , (30) 



Fi{z) = zW 



l + z 



Wia) - W 



1 



(31) 



Proof. We start with the function ^o(^)- To apply Theorem ll.4l we note that Fq{z) > for real z > 
(a e (0, e]) and Fq{z) is a holoniorphic function in the cut plane C\(— oo, —1 — l/(ae)] (cf. the branch 
cut B). For convenience, we define a function V{z) = 3Fo(z), then it remains to show that V{z) < in 
the upper half-plane. Since V{z) is a harmonic function in the domain > 0, it is subharmonic there. 
Thus we can apply either the maximum principle for harmonic functions in the form of [14, Corollary 
1.10] or the maximum principle for subharmonic functions (TSl p. 19-20]. In both cases, to get the desired 
result it is sufficient to ascertain that the superior limit of V{z) at all boundary points including infinity 
is less than or equal to f6|. In other words, V{z) < for 5z > if (cf. [HI p. 27]) 

lim V{z) < (3z > 0) 

and 

lim sup V{x + iy) <0 for aU x eR . (32) 



Since Fo{z) ^ 1/lnz for large z owing to ©, V{z) — > as |2;| — J> oo and the first condition is satisfied. 
To verify the second condition we note that V{z) is continuous (from above) on the real line z = x £ W. 
Then after introducing variable t = a(l + x) the inequality p2p becomes 




[W{t) - W{a) 



W{t) 



< 



or substituting W{t) — u + iv and using pUj) 



—V 



( 



V 



,2 



[(6 + wcot w)2 + w^]^ 



• 2 

sm V 



(33) 



where6 = M^(a) e (0,1]. 

We have = for t > -1/e and < t; < tt for -oo < t < -1/e (see ([3])). Since < 6 < 1 and 
z;^/sin^w > 1 for v e (0, tt), the inequality ([55)1 holds for all real t. Thus Fo{z) is a Sticltjes function. 
The theorem for the function Fi{z) follows from the relation 



because in terms of the conditions of Theorem ll.4l the transformation in the right hand-side of (IMl) retains 
the properties of Fq{z). In particular, '^Fi{z) < for 3z > because, first, and S>(— z/(l + z)) are 
of the opposite signs and secondly, SsFo(z) > for SJz < which follows from Fq{z) = Fq{z) due to near 
conjugate symmetry and the established above non-positivity of 3Fo(z) in the upper half-plane. Thus 
Fi{z) is also a Stieltjes function. 

Remark 2.6. We make a note about a behavior of functions pO|) and ()31|) for large and small z. 
Specifically, using 1^ and ([5]) one can obtain respectively Fq{z) — > and Fi(z) a/W'^{a) as z —> oo. 
Using ([T^ we find Fq^i ^ c/z as z 0, where c — {1 + W{a))^ /W{a). By the way, owing to property 
(0) given in theorem \1.2\ the differences i^o.i ~ c/z are Stieltjes functions as well. 



The principal branch of the Lambert W function itself is not a Stieltjes function in the sense of Definition 
11.11 This can be shown in different ways. For example, one can apply Theorem 11.31 to W{z) to see that 
the second condition (|23l) fails. Indeed, when z = is we have by ((T2)) and 



The same conclusion can be reached using Theorem II .41 because (jlj contradicts (|24|) . Finally, W is not 
a Stieltjes function because it is not an anti-Herglotz function (cf. Remark 1 1. 5 1) . 

Note, however, that W function can be regarded as a Stieltjes function in the sense of a definition 
given in [18] and [19] or used in [20] and different from pO| by the factor z in the right hand-side. W 
function can also be considered as a generalized Stieltjes transform by the definition in [2V (which is 
different from that of the generalized Stieltjes transform defined in [251 P- 30] and studied, for example, 
in [33] and [H]). Finally, in [3], the terms Stieltjes function and Stieltjes representation are not treated 
as equivalent (compare definitions f^, p. 11] and P", p. 55]). By these definitions W{z) has a Stieltjes 
representation (which is the result of multiplication of the representation ((77)) by z) though it is not a 
Stieltjes function. 

3 Particular Stieltjes integrals 

The Stieltjes representation for W/z given in ([27]) and (|28p itself contains W, which can be regarded as 
self-referential. Here we give representations containing only elementary functions. 

Theorem 3.1. The following representation of function W{z)/z holds f/^5| / 




(34) 



2.3 Is a Stieltjes function? 



|siy(is)| = s\/u'^ + = sec^(w)e 



V tan V 



— > cx) as V ^ -K 12. 




(35) 
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Proof. We start with 1^5)) and, noting change to the variable v = '^W(t). The integral becomes 

Wiz) _ 1 r V dv 



TT Jo t{z - t) V'{t) 



(36) 



where the variables t and v are related by The derivative v'{t) is conveniently found by taking the 

imaginary part of p^ . and using pHI) . 

^ " t[f2 + (l+u)2] = i[i;2 + (1 - WC0tw)2] ■ ^^''^ 

Remark 3.2. Since the integrand in (I35|) is an even function (with respect to v), the integral admits the 
symmetric form 

Wiz) 1 r w2 + (1 - t;cotl>)2 , , , 

— ' -dv (|argz| < tt) . 



z 2tt z + vcsc{v)e~'"'^°^^ 

This integral has a C°° periodic extension and thus the midpoint rule is spectrally convergent for its 
quadrature (see e.g. \46^ )- 

We now pay attention to the statement (d) of Theorem 12.41 with c = 1 which means by ([T4| that the 
derivative W'{x) G S. We derive an integral representation of W'{z) in the complex z-plane in the next 
theorem. 

Theorem 3.3. The derivative ofW function has the following Stieltjes-integral representation 

w'(^) = - r — — f\ -.cot. (larg^i < ^) • (38) 

Proof. We take the formula with a = due to p^ . After changing t by —t and introducing the 
function p(t) ~ —fi(—t) we obtain 

Wiz) = r ^ . (39) 



— oo 



The function p{t) can be determined using the Stieltjes-Perron inversion formula |25|, p. 591] 

p(t) = lim 3 / W'(t + is)dt 

TT s^Q+ 

for all continuity points t. Since p(t) is defined to arbitrary constant, after integrating one can set 

p{t) = -- lim 3Ty(i + is) ^ ~-^W{t) , (40) 

TT s-!-0+ TT 

where the limit uses the continuity from above of W on its branch cut. The same result can be obtained 
using one of Sokhotskyi's formulas 26,, p. 138]. 

Since p{t) increases from —1 to as i changes from — oo to — 1/e and vanishes on the remainder of the 
real line (cf. ([3])), the domain of integration in p9p is defined by — oo < t < —1/e. Owing to the relation 

^l{t) = -p{-t) , (41) 

the function p{t) in psp can be regarded as a positive measure. Thus takes the form 

TT J_oo t- Z dt 

This formula can also be found by considerations similar to those used in the proof of Theorem 12.11 By 
the change of variables v = ^W{t) in the integral (|^^ (see the proof of Theorem 13 . 1 1 and relation ([TT]) ) 
we obtain 



sm.^,,„, I d,^^, (43) 



Corollary 3.4. 



Proof. The integral ([55)1 can be written as 



1 r dv 



z-t 



(44) 



where t is defined by ([TT|) and the left-hand side is obtained by differentiation of the series (jH) that 
is convergent for \z\ < 1/e. Since |t| > 1/e and therefore \z\ < \t\, we can expand {z — t)~^ in the 
non- negative powers of z. Equating the coefficients of the same power of z in (|44p we obtain an equality 



-1) 



which after substituting pT|) results in (|43p . 



n-l ' 



dv 

7" 



It is obvious that if the integral is known then going back from it to (jH]) we find l|38p . The 
integral P5|) was conjectured by Nuttall for v > [27 ; Bouwkamp found a more general integral |5S], 
for which Nuttall's conjecture is a special case, using a representation of 7ri/''/r(i/ + 1) via a Hankel-type 
integral. Thus a representation of the derivative of W function in the form of the Stieltjes integral (l38l) 
allows one to compute the integral (|43p and conversely, starting with the integral of Nuttall-Bouwkamp 
one can obtain formula (|38|) in a way completely different from that used in the proof of Theorem [331 It 
is interesting to note that the connection between (j43| and Lambert W was noted by W.E. Hornor and 
C.C. Rousseau before W was named (see editorial remarks in 27 ). 

Coming back to the results of Theorem 12.41 we consider the assertion (a) with c = 1 and assertion 
(e) by which 1/(1 + W{z)) e S and 1/W{z) — 1/z G S. We can derive integral representations of these 
functions in the same manner as it was done for W'{z) in the proof of Theorem 13.31 The result is in the 
following theorem. 

Theorem 3.5. The following Stieltjes-integral representations hold 

1 I r dv 



1 + W{z) 



1 

TT 



'0 



1 



ze 



v cot V 



ainv/v 



argz| < tt) , 



(45) 



1 1 1 r -y^ + (1 -^;cot^>)2 , , ^ 

T77T^ = - H — / r-n ^-^ :—:dv argz < tt) . 

W{z) z TT 7o vcsc{v){vcsc{v) + ze""""^-") *= i ^ 



(46) 



Corollary 3.6. 



Wiz) = In 



^ ^ z 



v"^ + {1 — V cot w)^ 



TT Jq u csc(w) (f csc(u) 
Proof. By substituting gS]) in W{z) = \n{z/W{z)) (cf. ©). 



dv 



(47) 



4 Completely monotonic functions 

We denote by CA^ the set of all completely monotonic functions, which are defined as follows [12] . 

Definition 4.1. A function f : (0, oo) — > R is called a completely monotonic function if f has derivatives 
of all orders and satisfies (— l)"/*^"-'(a;) > for x > 0, n — 0,1,2, ... 

The set of Stieltjes functions is contained in the set of completely monotonic functions, and thus all of 
the functions listed in Theorem l2.4l are completely monotone. The set CA4 is a convex cone containing the 
positive constant functions; a product of completely monotonic functions is again completely monotone 
[SJ p. 61]. By Bernstein's theorem [21 Theorem 9.3], a function / e CAi if and only if it is of the form 

f{x) = / e--^dv{0 {x > 0), (48) 

where v is an uniquely determined positive measure on [0,cx)). Completely monotonic functions are in 
turn connected with the set of Bernstein functions. 



Q 



Definition 4.2. /5, Definition 5.1] A function f : (0, oo) — [0, oo) is called a Bernstein function if it is 
C°° and f is completely monotonic. 

We denote the set of Bernstein functions by B. Since W' E S C CAi, is a Bernstein function. A 
Bernstein function f{x) admits the Levy-Khintchine representation 



where a,b > and v is a, positive measure on (0, oo) satisfying ^(1 + S,)^^di^{S^) < oo. It is called the 
Levy measure. The equation (|49)) is obtained by integrating (|48)) written for /' 3 . 

An important relation between the classes S and B is given by the assertion |3i Theorem 5.4] 



Combining this with the function composition result |3, Corollary 5.3] that / G CM and g E B implies 
f o g E B, we obtain the following lemma 

Lemma 4.3. If f E CM and g eS\ {0} then /(l/.g) E CM. 

This lemma extends the list of completely monotonic functions containing W. 

Theorem 4.4. The following functions are completely monotonic 

(a) x^Wix) (x>0, X< -I). 

(b) x^W"'{xl^)[l + W{x'^)y (x > 0, a,7 > 0, -1 < /3 < 0, A < 0). 

(c) x^W"{x-l^)[l + W{x-P)y (x > 0, a,7 < 0, -1 < /3 < 0, A < 0). 

(d) l-a;^"^Tl^"T(a;'^)[l + iy(a;'^)]^^i (a; > 0, < a < 1, -1 < ;3 < 0, < 7 < 1). 

Proof, (a) Since W{x)/x E S C CM and x" E CM for a < 0, the function x^W{x) (A < -1) is a 
product of two completely monotonic functions and the statement (a) follows. 

(b) Take function fa{x) = x^"' E CM (x > 0, a > 0) and functions g{x) = l/W{x^) and h{x) = 
1/(1 + W{xP)) where -1 < /3 < 0. Since 1/g E S and 1/h G 5 by Theorem [23] (f) with c = and 
c = 1 respectively, by Lemma 1431 we have fa{g{x)) = g~°'{x) E CM and fj{h{x)) — h~'^{x) E CM 
(7 > 0). Substituting functions g{x) and h{x) in the power functions and taking a product of 
obtained completely monotonic functions with x^ E CM (x > 0, A < 0), the statement (b) follows. 

(c) Consider function f\{x) ^ x^ E CM (x > 0, A < 0) and functions g{x) — W{x~^) and h{x) = 
1 + W{x-I^) where -1 < /? < 0. Since l/g E S and l/h G 5 by Theorem [23] (g) with c = and 
c—1 respectively, by Lemma 1431 we have fa{g{x)) — g°'{x) E CM and fj{h{x)) == h'^{x) E CM 
for a < and 7 < 0. Substituting functions g{x) and h{x) and taking a product of obtained 
functions with f\{x), the statement (c) follows. 

(d) By Theorem [23] (h) and the assertion for a: > 0, < a < 1, -1 < /3 < 0, < 7 < 1 we 
have /(x) = g°''^{x)[l + W{x'^)]'-'-'^ E B, where g{x) = x-'^W{x'^). In addition, the function f{x) is 
bounded, particularly, < f{x) < 1 because < [1 + W{x^)]''^'^ < 1 and < g{x) < 1 (the latter 
follows from the fact that g{x) goes to and 1 as a: tends to and 00 respectively and g'{x) > 0, 
which can be estabhshed using ([S]), (0 and (HU)). Then by ^3, Remark 5.5] the assertion (d) follows. 

We considered only sufficient conditions for a function to be a completely monotonic. To find the 
necessary and sufficient conditions is a much more complicated problem so that in some cases it requires 
(at least as the first step) using the methods of experimental mathematics [29] . 





(49) 



gES\{0}^l/gEB. 



(50) 




in 



5 Complete Bernstein functions 



A very important subclass in B is the class of complete Bernstein functions denoted by CB. 

Definition 5.1. 14' Definition 6.1] A Bernstein function f is called a complete Bernstein function if the 
Levy measure in (|49p is such that dv{t)/dt is a completely monotonic function. 

We point out four connections between classes CB and S used in this paper (for additional relations 
between these classes see [H Chapter 7]). By Proposition 7.7 in [4], 

f eS^ fiQ)- f{x)eCB , (51) 

where the limit of f{x) &i x — Q (from the right) is assumed to be finite. Also if / is bounded and / G CB, 
there exists a bounded g £ S with lim2._j.00 gi^) = such that 

f{x)^f{0)+g{0)-g{x). (52) 

In addition, [H Theorem 7.3] and [4] Theorem 6.2(i),(ii)] cstabhsh 

f eCB^l/f eS\{0} , (53) 
f €CB^ f{x)/x e S . (54) 

Now we go back to the properties of the set S listed in section 11.21 to prove the last three properties 
therein. Let / e 5\{0}. 

(x) Apply sequentially (vii), (i), to obtain /" G 5 (0 < a < 1) ^ /"(O) - /"(x) e CB ^ 

g{x) = [r(o) - r{x)]-^ G 5 =^ i/.9(iA) = r(o) - r (1/x) &s- 

(xi) Apply sequentially dST]), (El), (ii), to obtain /(O) - f{x) € CB ^ g{x) = [/(O) - f{x)\ ^ G 5 ^ 
l/{xg{x)) = (/(O) - f{x))lx G 5 ^ (1 - f(x)lf{Q))/x G 5; 

(xii) By (vii), /" G 5 (0 < a < 1). Suppose that Vmvx^o f{x) = b < 00 and lim^^-^oo f{x) = c where 
< c < 00. Then < < c"°, i.e. is bounded. In addition, G CB by 1^. Therefore 
the statement (|52p can be applied, i.e. there exists a bounded function g G S,\imx^oo g{x) ~ such 
that we can write g{x) = g{0) + — f~°'{x). Taking the last equation in the limit x — > 00 we obtain 
g{0) + b~°' = c~", hence g = c~" — and the assertion follows. 

In closing this section we note that the statement ([55)1 with 1/W G S (by Theorem l2.4f a) with c — 0) 
immediately results in W ^ CB. Being a complete Bernstein fucntion W has an integral representation 
that is the result of multiplication of (pSj) by z [i', Remark 6.4], which reflects the relation (|54p . In 
addition, the complete Bernstein functions are closely connected to the Pick functions considered in 
Section [71 



6 Bernstein representations 

Not only does W G CB as shown, it also belongs to another subset of Bernstein functions. 

Definition 6.1. 14-' Definition 8.1] A Bernstein function f is called a Thorin- Bernstein function if the 
Levy measure in (|49p is such that tdv{t)/dt is a completely monotonic function. 

To find out whether ly is a Thorin-Bernstein function we apply Theorem 8.2 in [3], which establishes 
five equivalent assertions (i)-(v) which we refer to below. In particular, in accordance with assertions 
(i) and (ii), W{x) is a Thorin-Bernstein function because W{x) maps (0,oo) to itself, W{Q) = and 
W'{x) G iS. As a Thorin-Bernstein function, W{x) admits the following integral representation. 

Theorem 6.2. The principal branch of the W function can be represented as the integral 

VK(z) = ^ In (^l + z^e"™*"^ (|argz|<7r). (55) 

Proof. It follows from the proof of Theorem 8.2 in [4j and its assertion (v) that W{x) being a Thorin- 
Bernstein function, it can also be written in the form 

W{x) ^a + bx+ r -^^dt , (56) 
Jo x + t t 
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where due to W{Q) = and © we have a = and 6 = 0. Comparing and ^7^, (gll) we find 
T{t) — '^W{—t)/TT. Integration by parts gives 

W{x) = - ^ / In (l - I) j^^W{t)dt , (57) 

where we changed t — > —t. Changing the variables v — ^W{t) in (I57p (see and taking a holomorphic 
extension of the result to the cut z-plane C\(— cxd, 0] satisfying near conjugate symmetry, we obtain (|55p. 

Remark 6.3. Differentiating the representation (|55p for W{z) gives formula (j38p /or W'{z). 

Remark 6.4. The representation (|57p iwas obtained in ISOl as a dispersion relation for the principal 
branch of W function using Cauchy 's integral formula in a manner similar to the method applied for the 
proof of Theorem [ 



As a Bernstein function, W can be written in the form , with a = 6 = as in ([551) , allowing us to 
establish one more representation of W . 

Theorem 6.5. For the principal branch of the W function the following formula holds 

Wiz) = / ipiOd^ (3?z > 0), (58) 



^ - f exp {-Cv csc(w)e-'' ™* dv. 
^ Jo 



where 



Proof. We consider the Stieltjes integral (|T8l) for the derivative 

dn{9) 



(59) 



W'{x) -- 

and use the representation [x + 9)^^ = e^'^^+^-'^d^ to write it in the form 



oo f poo 



W'{x)=J |y e-^'^dfi{e)\e-''^d^ . (60) 



Comparing and the result of differentiating P^]) wc find the relation between measures /x and i/ [21] 

e-^^diiiO) . 



dv _l 

Using the relation (|4T|) and formula (|40)). after changing the variables v = '^W{—9) (see (fTTt ) we obtain 

dv = ^d^ , (61) 

where (p{^) is defined by ([5^. We collect the intermediate results and take a holomorphic continuation 
of (|3H]) to the right half-plane > 0, where the integral ([55| is convergent, in accordance with near 
conjugate symmetry (cf. Proposition 3.5 in [3]). 

Note that by (|59l) function (y5(^) G as should be, because W is still a Thorin-Bernstein function (cf. 
Definition 16.11) 



7 Pick representations 

Definition 7.1. ^ Definition 4-1] A function f{z) is called a Pick function (or Nevanlinna function) 
if it is holomorphic in the upper half-plane > and S/ > there. 
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A Pick function f{z) admits an integral representation [3] Theorem 4.4] 

f°° 1 + tz 

f{z) ^ao + boz+ ——-da{t) {Qz > 0) , (62) 

where 

ao = 6o = lim ^ , (63) 

and a positive measure a satisfies 

lim - I '^f{t + is)(p{t)dt = I Lp(t)dcr{t) (64) 

for all continuous functions (/? : M ^ R with compact support. The formula ()62p with the integral written 
in terms of a measure da{t) — 7r(l + t'^)~^da(t) is called a Nevanlinna formula [32^ p. 100]. 

Since T4/^(z) is a holomorphic function in the upper half-plane > with the property W{z) is 
a Pick function. It also follows from the two facts that W G CB (see Section [5]) and that the complete 
Bernstein functions are exactly those Pick functions which are non-negative on the positive real line [31 
Theorem 6.7]. Thus W admits a representation ([62]) and in view of that the following theorem holds. 

Theorem 7.2. The principal branch ofW function can be represented in the form 

1 r 

W{z)^aa + - K{z,v)t{v)dv (|argz| < tt) , (65) 



where = ^W(i) = 0.3746990. 



{l + zt{v)){v^ + {l~vcotvf) 

K{z,v) = , 66 

(z-t(t7))(l + t-(w)) 

and t{v) is defined by pip . 

Proof. Apply formulas (jMll-dMl) to function f{z) = W{z) 

W(iv) 1 
ao = ^W{i), bo = lim — da{t) = -^W{t)dt. 

j/->oo iy TT 

Using dll), we see bo = 0. Since ^W{t) = for t > -1/e (cf. ©), we obtain 

W{z) ^ao + - f W(Orft (3z > 0) . (67) 

i-oo (i-z)(l+t^) 

By the change of variables v = '^W(t) in the integral (1571) (see PTjl ) we obtain formula that is also 
valid in the lower half-plane 3z < in accordance with near conjugate symmetry of W . 



^^=7oexp(-i/ K{z,v)t(v)dv\ (|argz| < tt) , (68) 



Corollary 7.3. 



where 70 = e-'^^t^) = 0.6874961... 

Proof. It immediately follows from (|65l) owing to the identity Ty(z)/z = e"'^^^). 

Now we take advantage of the fact that if function f E S then — / and 1/f are Pick functions [3]. 
Therefore, since W{x)/x S S, —W{x)/x and x/W{x) are Pick functions that admit a representation 
([62]). We can obtain a representation (l62|) for functions — iy(a;)/a:: and a;/iy(x) similar to the derivation 
of formula (I65|) . and the result is in the following theorem. 
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Theorem 7.4. For the ■principal branch of the W Junction the following formulas hold 

W(z) 1 f 

+ - K{z,v)dv (|argz|<^), (69) 

Z Jo 

^ m-- r K{z,v)e~^"'°'^dv (largz|<^), (70) 



W{z) TT _ 

where K{z,v) is defined by /3o ^^[W{i)/i] = QW{i) = 0.5764127.., r]Q = ^[i/W{i)] = 1.21953... 

The constants in ([65t and ((68)l - ((70|) obey the relations + ifio = W{i), 70 = e^"° — Po/ cos/3o, 

We add in one more integral representation associated with the Nevanlinna formula which follows 
from the result obtained by Cauer jSSJ . Specifically, based on the Riesz-Herglotz formula [351 p. 99] 
Cauer proved that if a real symmetric function f{z) with non-negative real part is holomorphic in the 
right z-half-plane, it can be represented as 

dh{r) 



f{z) = z 
where constant 6 > and 



(71) 
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h{r) = - lim SR / fix + iy)dy . (72) 

TT x^O Jq 

In fact, the formula ((7T|) follows from the Nevanlinna formula (or (|62p ) after changing the variable 
z —iz, which transforms the upper half-plane onto the right half-plane, and taking into account 

fi-z) = m- 

Theorem 7.5. The following representation of function W{z)/z holds 

W{z) 2 f^/^ [v^ + il + vtai-iv)']vsec{v)e^'-''\ ^ n^ 

— ' tSLUvdv [Uz > 0) . (73) 



z TT Jq -I- ^2 sec^(w)e^"* 

Proof. Since W function meets the above requirements, the formulas ([7T|) and (17^ can be applied with 
the result 

W{z) 2 f°° ^W{is) 



2 f°°^W{is)^ . 



where we set 6 = due to © and r = . 

Changing the variables defined by (|12l) - p^ we obtain 



W(z) 2 i;tanw ds , 

— ' dv . (74) 



z TT 7o + s'^{v) dv 
Similar to ([57| one can find 

dv V , , 

ds s{v)[v'^ + (1 + vta.nv)'^] ' ^ ' 

Substituting ([75)1 and into (fTl)) . the theorem follows. 

Remark 7.6. Comparison of the formula (|73p with the representation (j35p (taken in the right z-half- 
plane) shows that the integrand in the former contains z^ rather than z, which can he profitable in using 
the integral representations for numerical evaluation ofW{z) at large z. 



8 Poisson's integrals 

Theorem 8.1. The following two formulae of Poissot^ [34l hold for x G (— l/e,e) 

2 r cos|6'-a:e-™'^''cos(|6'H-sin6') , , , 

W(x) = - / — ^ 5 5 — 5^ — ^ cos d0 76) 

^' TT L l-2a;e-™''''cos(6i-t-sin6l -)-a;2e-2cos0 2 ^ ' 



^The second formula is explicitly given in |34l sec. 80, p. 501] in terms of the tree function T(x) (see I I82I I and I I83I I 
below) and proved using the Lagrange Inversion Theorem 1351 p. 133] and a series expansion of the logarithmic function 
— ln(l — e*^(/)) in powers of e'^ where the expansion coefficients 4>" /n are exactly the coefficients of the complex exponential 
Fourier series for the same function. On the other hand, today it is well known 1361 p. 143-145] that there is a tight 
connection between the classical Poisson Formula and the Cauchy Integral Formula. Our proof is based on the latter and 
thereby differs from that given in the original. 



1/1 



2 r sin 1 6* H-se™'''' sin (1 6* -sin 6, , , , 

Wix) = / — ^ At;^ , J „ sin U AO 77 

' n In l + 2a:e™''«cos 6'-sin6') +a;2e2™''e 2 v ^ 





Proof. We consider the defining equation ([T]) for given real z — x 

We^ = X (78) 

and interpret it as an equation with respect to W. Then we can write the equation in the form F{W) — 
where 

^^(C) = C - a:e-^ . (79) 

Let r be the positively-oriented circumference of the unit circle |C| — 1 in the complex C-plane and 
domain G be the interior of F. The function F{C) is holomorphic in G and by Rouche's theorem it has 
a single isolated zero there when \x\ < 1/e because in this case |— a;e~^| < |C| on F. Therefore, using 
Cauchy's integral formula with taking F for the integration contour we can write 

for |a;| < 1/e. 

Since F'(C) = 1 + xe^^ = 1 + ( by ([711) and we obtain 

W^i^ ^ t ,,1 . ,, de , (81) 

2ji J \ 2,g— cos 0— i(e+sin 9) ' ' 

where we set C — e*^,— tt < 9 < tt. Separating the real and imaginary parts of the integrand in (1811) 
we find that the former is an even function of 9 whereas the latter is an odd one. Thus, the integral of 
the imaginary part vanishes, as should be, and the integral of the real part gives double the value of the 
integral on [0,7r]. As a result, after some arrangements, we come to integral ([76]). 

If, instead of ([78| . we consider the equation defining the (Cayley) 'tree function' T{x) |37l p. 127-128] 

Te"^ = X (82) 

and introduce function H{Q — C,~xe'^ in a similar way as function (|79p then after analogous considerations 
and taking into account a relation 

W{x) = -r(-x) (83) 



in the final result we obtain formula (|77|) . 

Now we discuss the domain of validity of the integrals (|76|) and (|77p which is actually wider than the 
interval —1/e < x < 1/e arisen above in applying Rouche's theorem. It immediately follows from the 
fact that is a single valued function and therefore i^(C), the denominator in (|5(I1) . has a single zero in 
G for each such x that |C| < 1, i-e. for — 1/e < x < e. Since Rouche's theorem is a consequence of the 
argument principle (see e.g. [38]), it would be instructive to obtain this result using the latter. To do this, 
say for integral ([76]) . we apply the argument principle to function ([79| in case when a; > 0. It is easy to 
see that function rj = F{Q performs a conformal mapping of the strip {— oo < 3?C < cxd, — vr < 3C < tt}, 
containing the entire domain G, to the complex //-plane cut along two semi-infinite lines on which rj = 
5 ± ZTT,^ > 1 -I- Incc. We also cut the ?7-plane along the negative real axis to take |arg?7| < tt in the cut 
plane and consider an image of F which is defined by equations 

pcosip = cos 9 - xe^ ""'^ ^ cos(sin 9) , (84a) 
p sin ip = sin 9 + xe" ^ sin(sin 9) , (84b) 

where p = |ry| and ip = argry. 

The equations (|84|) are invariant under transformation 9 — ;> —9, ip ~p and describe a closed curve 
F that is symmetric with respect to the real axis in the ?7-plane. Suppose that while a variable point ^ 
moves along F once in the C-plane, the image point r] = F{C,) moves on F once in the ry-plane, making 
one cycle about the origin. Then the change in argument of rj is 27r and therefore, by the argument 
principle the function F{C,) has a single zero in G [38] p. 48]. For this it is necessary that two points on 
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r corresponding to 1^9 = and ip — -k are located on the real axis on the opposite sides of the origin, i.e. 
with positive p to be measured on the opposite rays. Substituting = tt in (j84l) gives pcomp — --l — xe 
and psinif = 0. It can be p > only when ip — t:; then p = 1 + xe is positive for any a; > 0. When 
9 = 0, we have p cos 95 = 1 — x/e and psimp = 0. Now ip = and p = l— x/e>0 when x < e. Thus 
for < a; < e the curve F encloses the origin. Since for these x the right-hand side of equation (|84bl) 
vanishes, i.e. ^rj = sequentially at = —n^O = and 9 = tt as 9 continuously changes from — tt to tt, 
the curve F is traversed once with exactly one cycle about the origin being made. This corresponds to 
the fact that the inverse of the mapping 77 = ^(C) is continuous in the domain bounded by the curve F 
and on F itself and hence F consists only of simple points '39', Theorem 2.22]. Thus, by the argument 
principle the function -F(C) has a single zero in G. Summarizing up the obtained results we conclude that 
the integral (|76l) is valid for x G (— 1/e, e). The integral ((77)) can be considered in a similar manner. 

Remark 8.2. The integral representations (|76p and (j77p can be immediately applied to the tree function 
using relation 



Remark 8.3. We can apply the above approach to the equation To eliminate a singularity at the 
origin we compose the integration contour of a small circle of radius, say r, and the unit circle, both 
centered at the origin and connected through the cut along the negative real axis. Then, making r go to 
zero, we find for < x < e 

9 /, • 3 /, 3 
cos — h fsm — (/ — cos — wma; n 



where 



W{x) = ^{x) + I _ ^ , , ^„ , , 2 cos "-d9 

f-i 



l + 26isin6' + 02_2cos01na; + ln^x 2 



ibix) = / ; — —dt . 

lo 7r2 + (lna; + i-lnt)2 



9 Burniston-Siewert representations 

One of the analytic methods for solving transcendental equations is based on a canonical solution of a 
suitably posed Riemann-Hilbert boundary- value problem [26, p. 183-193]. This method was found and 
developed by Burniston and Siewert [S], its versions, variations and applications were also considered by 
other authors. By the method, a solution of a transcendental equation is represented as a closed-form 
integral formula that can be regarded as an integral representation of the unknown variable. Below we 
consider such integrals for the W function which are based on the results of application of the Burniston- 
Siewert method to solving equation (|78|) obtained in paper [4^ and the classical work [42] . 
We start with two formulas derived in [3^ and apply them to function (fZU]) 

..W^-F(0)»p{-^/™<l.c}, ,85, 



where the integration contour F is the unit circle \(\ = 1 and x G (— l/e,e). Since F(0) = —x and 
W{x)/x = e-'^(^), formula §B is simplified 

Wix) = ± I M£iOZO . (S7) 



27rj Jy C 

We set C = e*^ -tt <9 <ti. Then, as F(C)/C = F(e^'')e-*^ = R{9) + il{9), where 

R{9)= l-xe-""''^ cos{9 + &ui9), 
I{9)= a;e-™'^^sin(e' + sin6l), 

and dC/C — id9, the integral ([57]) is reduced to 

W{x) = ^ f In [R^[9) + /2(6i)) d9 . 







Similarly, the integral ([55]) can be represented in the form 
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W{x)^—l {2 aictanil (9) /R{0)) sin 0~ In {R^ (9) +1^ (9)) cos 9} d9 , (89) 

where we have taken into account that aig{R{9) + il{9)) = aTctan{I (0) / R{9)) as R{9) > for < < tt 
and — 1/e < x < e. We note that the integral (|88l) has a simpler form than ([M)) . Integrals similar to the 
above with using a function $(C) = C^^ — a; in our notations instead of ([7^ in formulas (|55]) and (1551) 
(without simplification (|57|) are given in (3D]. 

Thus the integrals (j88l) and (|89)) representing the principal branch of the Lambert function are 
valid in the domain that contains interval (— 1/e, 0). However, there is one more branch that is also a 
real- valued function on this interval, this is the branch —1 with the range (— oo, — 1) (recall Wq > — 1 and 
Wo(— 1/e) = W_i(— 1/e) = —1) [1 . A representation of this branch can be obtained on the basis of a 
simple interpretation of formula (j86|) given in |41) 
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W.,ix) = l-2c-—J^lni^-^jdC, (90) 

where the circle C is defined by equation |C + c| = c — 1 with arbitrary constant c > 1 and — 1/e < a; < 
— (2c — l)e^~^'^. Transformations of ([TO]) leading to a definite integral are similar to those used above to 
obtain the integrals ([55)1 and ([M)) and skipped here together with a bulky result. 

We return to the principal branch and use the result in |33J formula(13)] to write [33] 

»' w - - (i„ . - 1) f >„ ( i:::;:;:;:;: ) m 

or 

f 1 f°° arg(ln z + t — Int + iir) ,1 , , 

W^(z) = l + (lnz-l)exp|--y^ ^-p^ '-dtj, (92) 

where z ^ [—1/e, 0]. In case of real z = a; > 1/e, when the expression Inz + 1 — int is real and positive 
(for t e (0, oo)), the formula (j92l) is simplified and reduced to 

W{x) = 1 + {inx -l)cxpl-^ [ arctan ( ^ ^ /[ — T-t) "TT-tI (93) 
or, after integrating by parts 



TT Jo \in X + t ~ Int J 1 + t 



We emphasize that the domain x > 1/e of validity of the formulae and ([M)) is different from that of 
(pl and (p)). 

For the case a; G (—1/e, 0), we refer the reader to [311 formula (32)] where the principal branch Wq and 
the branch W-i are represented in the form of a combination of two expressions similar to the right-hand 
side of (|M]) . 

Remark 9.1. We can regard the integral in the formula (|9ip as an improper integral depending on 
a parameter p ~ inz and consider it in the limit p — > oo (when z ^ oo). Since the integrand is 
a continuous function of two variables t and p in the domain under consideration and the integral is 
uniformly convergent with respect to p, we can take the limit under the integral sign and find that the 
integral vanishes as the integrand goes to zero. Then the formula (PT]l reproduces the asymptotic result 

Finally we note that by use of elementary complex analysis in |44| there is obtained a common closed 
form representation for all the branches Wk{z) in the complex z-plane through simple quadratures. 
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10 Concluding remarks 

We derived various integral representations of the principal branch of the Lambert W function using 
different approaches. The most part of them is associated with functions containing W belonging to 
various classes of functions admitting certain integral representations. Among other classes we considered 
in detail the classes of Stieltjes functions and complete monotonic functions and by the example of 
functions based on W in fact demonstrated different ways to establish belonging of a function to these 
classes. 

Besides their own importance the derived integral representations have some applications. One of 
them has been mentioned in connection with finding Nuttall-Bouwkamp integral (I43p . Other definite 
integrals appear in taking the obtained integrals with a particular value of z. For example, integrals (j35p . 
(|i5)) . (gni), dZSl) taken at z = e yield respectively 



Another advantage that can be taken of the obtained results is based on a comparison between different 
representations of the same function. This reveals equivalent forms of the involved integrals. In addition, 
since some of the integrals are simpler than others, such equations can be regarded as a simplification 
of the latter. For example, equating integrals (|69p and (|55]) shows that the former can be simplified and 
reduced to the latter. 

At last we mention that the Pick representations (|65p . ((68|) . (l69t . and ((70)) can be considered as 
integrals expressing properties of the kernel K{z, v) defined by (j66l) . 
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